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ABSTRACT 

A X-ray background synthesis model is used to calculate the contribution of 
Active Galactic Nuclei (AGNs) to the 1.4 GHz number counts between 100 nJy 
and 10 mJy. The number counts are broken down into contributions from radio- 
quiet and radio-loud AGNs, obscured and unobscured AGNs, and for different 
ranges in redshift and 2-10 keV X-ray luminosity, Lx- Compton-thick AGNs are 
included, but only to the level required to fit the peak of the X-ray background. 
The predicted radio counts show that the /iJy AGN population will be dominated 
by obscured, radio-quiet Seyfert galaxies with logLx < 43, and spanning < 
2; < 3. However, depending on the exact relationship between the radio and X- 
ray luminosities in radio-quiet AGNs, additional radio flux due to star-formation 
within AGN host galaxies may be necessary in order to match the observed 
AGN counts at a flux density of ~ 50 yuJy. The star-formation rates (SFR) 
required are modest, only ~ 3 Mq yr~^, assuming a constant rate with z and 
Lx- A more observationally and theoretically motivated relationship, where the 
SFR cx (1 + z)^'™(logLx — 40)^-^, will also account for the observed counts. The 
fiJj AGN population will provide a very clean sample to trace the accretion and 
galactic star-formation histories of Seyfert galaxies over a significant fraction of 
cosmic time. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: Seyfert — 
galaxies: starburst — radio continuum: galaxies 



1. Introduction 



The identification of p o werful radio sources with quasi-st ellar objects in the 1960s 
( iMatthews fc Sandagd Il963l : [Hazard et al.l Il963l : ISchmidtl Il963l ) ultimately led to the re- 
alization that Active Galactic Nuclei (AGNs) were the locations of accreting supermassive 
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black holes (e.g., ISalpeteil Eoei : iLvnden-Belll Il969h . Since that time, radio observat 



ions 



have continued to be an important probe of the AGN phenomenon, especially in the inves- 



tigation of relativistic jets emanating fr om the central engine (e.g., iBridle &: Perleyl Il984j : 
Kellermann fc Owen! Il988l : lO'Deal |2002| ). The next decade will see a major leap forward 
of the capabilities of ce ntimeter wave radio instrumentation with the development of the 
Expanded- VLA (EVLA; Napier 2006 ) and, ultimately, the Square Kilometer Arra}0. The 
increase in sensitivity afforded by these new instruments will probe the radio population 
down to flux densities less than a /xJy, opening a new window into the study of AGN evolu- 
tion and the growth of supermassive black holes. 

At flux densities S > 1 mJy, the radio source population is dominated by radio-loud 
AGNs, and it has long been known t hat the numb e r cou nt distribution indicates strong 
cosmic evolution (e.g., iLongair Il966l : ICondonI Il984j . Il989h. However, the radio number 



counts significantly flatten at fluxes below 1 mJy (e.g., 



Hopkins et al.lll998l : lRichardsll2000l : iSeymour et al.ll2004l : iSimpson et al. 



Condon 



1984: 



Windhorst et al 



2006 



1985 



Kellermann et al 



20081 ) indicating either a sudden change in evolutionary properties, or, more likely, the ap 



pearance of a new population of source s such as emission from star-forming galaxies (e.g., 
Windhorst et al.lll995l : iBenn et al.lll993l : lRichardsll2000l : ISeymour et al.ll2004l ). Strong radio 
emission from star-forming regions is expected as supernova remnants will accelerate cosmic 
rays which will radiate synchrotron emission in local magnetic fields. At low redshifts, the 
radio and far-infrared powers of star-forming galaxies are tightly c orrelated , allowing for the 



radio luminosity to act as a star-formation rate (SFR) indicator ( CondonI 1992 : Yuii et al. 



2001 


; iBelll 


2003) 


Ibar et al. 


20081; 



Belli 2003 ) . This correlation seems to hold at higher redshifts ( Beswick et al. 



2008 



In recent years, severa l surveys have cont i nued to probe the radio source population 
to ever-smaller values ( e.g. jHopkins et al.ll2003l : ISimpson et al.ll2006l : ISchinnerer et al.l 120071 : 
Kellermann et al.ll2008l ). with the deepest VLA surveys reaching rms noise levels of only a few 
yuJy per beam. Matching the radio sources to complementary optical. X-ray and infrared 
surveys has allowed the properties of raany of the faintest r adio sources to be estimated 
dSevmour et"aDl2008l : ISmolcic et allbood : IPadovani et"aDl2009h . The results of these studies, 
while still with large error-bars, indicate that radio-quiet AGNs rn ay make up to one-third 
of the radio source population below 100 /xJy (IPadovani et al.ll2009l ). As ~90% of AGNs are 
radio-quiet, these results imply that the deep radio surveys are now uncovering the dominant 
population of active galaxies. 

The last decade has provided a wealth of new data tracing the evolution of supermassive 



^http : //www . skatelescope . org 
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black holes and their impact on galaxy formation fe.g.. iKauffmann et al.ll2003l: iGrogin et al. 



20051: iPierce et al 



2007: Silverman et al. 2008a: Gabor et al. 



200S 



Schawinski et al 



2009 



Hickox et al.ll2009l ). but it has proved difficult to simultaneously study the AGN and host 
galaxy. Galaxies are much less luminous than AGNs at X-ray energies, so deep X-ray ob- 
servations are very efficient at fir iding all but the most heavily obscured AG Ns over a wide 
range of redshift and luminosity (IMushotzkyl |2004J : iBrandt &: Hasingerll2005l ). These surveys 
have showed that the vast majority of AGNs in the Universe are absorbed, and that AGNs 



evolv e in the same form of 'cosmic downsizing' as the galaxy population (e.g., lUeda et al. 
20031 : iBarger et al.l 120051). However, aside from indirect inferences using the evolution of 



the X-ray absorber (IBallantyne et al.ll2006l ). it is very challenging to glean much informa- 
tion about the affect of the AGN on the host galaxy through X-ray observations. Optical, 
ultraviolet, and infrared wavelengths allow study of both the host galaxy and AGN, but 
the extinction and re-emission properties of dust in the host galaxy severely complicate the 



Moran et al. 


2002: 


Donley et al. 


2008) 



nearly uncorrupted view of both processes. All AGNs produce some form of non-thermal 
nuclear rad fo emission which, although weak, can be identified by its high bri ghtness temper- 



ature (e.g.. iHutchings fc Nefflll992l : iBlundell et al.lll996l : iKukula et al.lll998l ) and is immune 



to attenuation by any circumnuclear torus. Radio-emission from the host galaxy will provide 
a similarly unabsorbed view into ongoing star-formation. Thus, combining studies of the /iJy 
radio-quiet AGN population with other multiwavelength surveys will allow a relatively clean 
look into the combined problem of AGN and galaxy evolution. 

In this paper, we make use of a X-ray background synthesis model to predict the number 
counts of radio-quiet AGNs down to a flux density of 100 nJy. The calculations break down 
the contribution of AGNs from different redshifts and nuclear luminosities. In addition, 
we add in radio-emission from various levels of star-formation in the AGN host galaxies 
to explore the impact on the predicted number counts. The next section describes the 
calculations necessary to move from a X-ray background synthesis model to predicting the 
radio number counts. Section [3] then presents and describes the results, including the effects 
of adding in star-formation. The results are discussed and summarized in Section 4. The 
following A- dominat ed cosmology is ass umed in this paper: Hq = 70 km s~^ Mpc~^, Qj>^ = 
0.7, and Qm = 0.3 (ISpergel et al.ll2003l ). The flux density of a radio source has a spectral 
index a, so that S* oc z/". 
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Calculations 



Previously, IJarvis &: Rawlinga (l2004j ) and IWilman et al.l (|2008l) est i mated the radio- 
quiet AGN contribution to the radio counts by employing the lUeda et al.l (120031 ) hard X-ray 
luminosity function, multiplying by a consta nt to correct for Compt on-thick AGNs, and 
then employing a i?05'^T-based correlation (IBrinkmann et al.l |2000| ) to convert from X- 
ray to radio luminosity. In both these papers, they use separate radio-derived luminosity 
functions to account for radio-loud AGNs. Below, we present a new method to estimate the 
radio-quiet AGN number counts that will improve on this previous work in several ways. 
First, we base our calculations on a fit to the X-ray background. This allows the peak 
spectral intensity of the background at ~ 30 keV to be used as a constraint on the missing 
Compton-thick fraction, and has the added feature of allowing us to break down the number 
counts into contributions from different groups such as obscured AGNs, or ones in a specific 
redshift or luminosity range. Furthermore, we account for a radio-loud fraction of AGN. 
The X-ray luminosity function counts both radio-loud and radio-quiet AGNs, so the radio- 
loud fraction must be subtracted from the total before predicting the radio-quiet counts. 
Finally, we convert from X-ray to radio luminosity using a luminosity-dependent conversion 
factor derived from high-resolution radio and hard X-ray observations of the nuclear regions 
of nearby AGNs. While clearly there remain several uncertainties in the steps described 
below, the changes outlined above should result in a more accurate prediction of the AGN 
contribution to the deep radio counts. In these calculations, an AGN is assumed to be 
obscured if it is attenuated by a hydrogen column A^h > 10^^ cm~^. 



2.1. The X-ray Background Synthesis Model 



The hard X-ray background encodes within it the entire history of accretion onto su- 
permassive black holes, independent of whether or not an individual AGN is radio-loud or 
radio-quiet. A detailed description of the X-ray background synthesis model employed here 
is beyond the scope of the current paper, and, indeed, many aspects of the model such as the 
X-ray spectral shapes, and the distribution of column densities have no impact on the pre- 
dicted radio number counts. For this application the important ingredients of the model are 
the absorption-corrected hard X-ray luminosity function, and the obscured-to - unobs cured 
AGN ratio (i.e., the ratio of type 2 AGNs to type 1 AGNs). The lUeda etaP J2003h hard 
X-ray luminosity function is the basis for the calculations presented here. 

The AGN typ e 2 / type 1 ratio is obse r ved to significantly decrease with AGN luminosity 
(lUeda et al.l l2003l : iLa Franca et al.l l2005l : ISimpsonI l2005l : lAkylas et al.l l2006l : iTreister et al. 
20081 ). The explanation for this effect is not entirely understood, but is most likely related to 
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the increased o utward radiation pressure produced at high AGN luminosities ( lLawrencelll991 
Simpsonll2005l ). The best determination of the local {z ~ 0) dependence of the AGN type 
2 /typ e 1 ratio with lumi nosity was produced by the Swift /BAT detector in the 14-195 keV 
band (ITueller et al.ll2008l ). AGN selection at these energies is expected to be completely unbi- 
ased to objects with Compton-thin absorbers. Figured] plots the BAT AGN type 2 fractions, 
/2 as a fu nction of Lx= L2-W keV where the conversion logLx ~ 1.06 log Li4_i95 kcv — 3.08 
found by IWinter et al.l (120091 ) was used to convert the BAT luminosities to Lx- The plot 



shows a strong trend of an increasing fraction of obscured AGNs as the nuclear luminos- 
ity decreases. The lowest luminosity point indicates a much lower obscured fraction than 
expected by this trend, but lies at such low luminosity that it is likely contaminated by LIN- 
ERs and other low-luminosity AGNs that have significantly different nuclear en vironment s 
and evolutionary histories than rapidly accreting AGNs at higher luminosities (iHol l2008l ). 
In practice, the hard X-ray luminosity functions are integrated down to logLx = 41.5. A 
power-law of the form /2 cx L^^ was fit to the data with the normalization set by the type 
2/type 1 ratio at logLx = 41.5. The lowest luminosity datum was included in the fit, but 
its removal makes a negligible impact on the results. Assuming an obscured-to-unobscured 
ratio of 4:1 at logLx = 41.5, the best-fit (reduced = 0.97) was found for (3 = 10.5. The 
predicted rad i o num ber counts decrease by only 6% if the smooth step function shown by 
Tueller et al.l (120081 ) is used to describe the variation in /2. 



This description of the AGN type 2 fraction is based on observations of local AGNs, 
but there is now accumulating evidence that also increases with redshift as (1 + z)'^ 
dla Franca et al.lbood : [Ballantvne eTaPlgOod: Frreister fc Urrvlbood : Easingeil l2008h . This 
finding is still controversial ( Dwelly fc Pagel 20061 ). and the data is not yet at the point to 
determine if there is a luminosity dependence to 7. Here, we assume 7 = 0.4, a value 
consistent with the various estimates, and that the AGN type 2 fraction evolves out to 



z = 1, where, in analogy with the cosmic star- formation rate density (e.g., iHopkind |2004J : 
Hopkins fc Beaconj |2006| ). it f lattens to 7 = 0. If we instead choose 7 = 0.62 with the 
evolution continuing to z = 2 (IHasingerl |2008| ) . then the predicted radio counts increase by 
only 8%. 

The hard X-ray luminosity functions derived in the literature are based on observations 
at energies less than ~ 10 keV, and so are insensitive to Compton-thick AGNs, those with 
absorbing column densities iVn > 10^^ cm~^. Therefore, the unknown fraction of Compton- 
thick AGNs has to be included 'by hand' int o any X-ray backg round model in order to fit 
the peak of the spectrum at ~ 30 keV. For the lUeda et al.l (120031 ) X-ray luminosity function, 
it was found that a Compton-thick fraction equal to the fraction of Compton-thin type 2 
AGNs provided a very good fit to the peak of the X-ray background. That is, at every z 
and Lx, the Compton-thick AGN fraction was equal to f2- The implicit assumption is that 
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the Compton-thick population evolves in the same manner as the Compton-thin AGNs. 



2.2. Translating the X-rays to Radio 



The X-ray background model provides the number density of obscured and unobscured 
AGNs as a function of z and Lx- To calculate how these AGNs appear in the radio sky, 
the X-ray luminosity must first be translated into radio lumi nosities, which, for radio- quiet 
AGNs, is not necess arily straightforward. Previous work by 
Wilman et al.l (120081 ) made use of the correlation published by 



Jarvis fc Rawlings (120041) and 
Brinkmann et al. ( 200ol ) based 



on AGNs selected through cross-correlating the ROSAT All-Sky Survey and the FIRST 
1.4 GHz radio catalogue. However, it is not clear if this correlation is valid over a wide range 
of AGN luminosities. The ROSAT All-Sky S urvey was a shallow, soft X-ray survey and 
the sources selected by iBrinkmann et al.l (120001 ) are typically unobscured, luminous quasars. 
These objects lie at substantial redshifts {z > 0.1) and therefore its possible that the radio 
emission in the FIRST beam may be contaminated by star-formation in the host galaxies. 
As we are interested in the relationship between the core X-ray and radio luminosities, it 
would be best to use results based on observations of lower- luminosity radio-quiet AGNs 
that have isolated the nuclear radio emission. 



Therefore, we consider the results of iTerashima fc WilsonI (120031 ) and iPanessa et al. 



( 120071 ) who studied the X-ray and radio nuclea r luminosities of nearby AGNs with high reso- 
lution observations. The iPanessa et al.l (120071 ) sample was limited to low luminosity Seyfert 
galaxies and LINERs with logLx < 43, with many of the AGNs accreting at very low rates 
with respect to the Eddington limit. This is probl ematic because it is now known that radio 



loudness is more common at low accretion rates (ITerashima fc Wilson! l2003l : IPanessa et al. 



20071 : ISikora et al.ll2007l ). so a s ample of low luminos i ty AG N may be biased toward high 
radio luminosities. In contrast, ITerashima fc Wilson J2003h plot the AGN radio to X-ray 
luminosity ratio as a function of Lx over a wide range of X-ray luminosities. Although there 
is a scatter of ~ 0.5 dex, they find that it decreases with Lx until it reaches a minimum at 
logLx ~ 43-44 before increasing again. Defining Rx = logz/Lj^(5 GHz)/ log Lx we quantify 
this behavior as 



R 



X 



-0.671ogLx + 23.67 41.5 < logLx < 43 

-5 43 <logLx<44 

logLx - 49 44 < logLx < 45 ^ ^ 

-4 logLx > 45 

This is the relationship used to convert from the X-ray to the rest-frame 5 GHz radio lumi- 
nosity. To calculate the observed flux at 1.4 GHz, a spectral index of a = —0.7 is assumed 
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( Kukula et al.lll998l ). In Sect. 13.21 we will show how the results change if the iPanessa et al. 



( 120071 ) relationship is used for all luminosities. 



The X-ray to radio conversion for the radio-loud population is also derived from lTerashima &: Wilson 



( 120031 ): z/L^(5 GHz) 10^^ Lx- This proportionality factor has a scatter of ~ 1 dex, but, as 
is seen below, the radio-loud population does not impact our conclusions, as they comprise 
a negligible component of the radio population at /iJy levels. A spectral index of a = —0.8 
is assumed to calculate the observed 1.4 GHz fluxes of radio- loud AGNs. 

The small fraction of radio-loud AGNs has to be removed from the entire X-ray selected 
population to concentrate on the properties of the dominant radio-quiet AGNs below 100 /xJy. 
As mentioned above, radio-loud AGNs are most common at high luminosity (the radio-loud 
quasars) and at very low nuclear luminosities, but these latter objects are not important to 
the X-ray background and are omitted here. At every z, the radio-loud fraction is assumed 
to rise exponentially from logLx = 41.5 to 46: /rl = cexp(logLx — 40) + b. A decent 
fit to the radio number counts between 1 and 10 mJy was found assuming /rl = 0.0175 at 
logLx = 41.5 and /rl = 0.10 at logLx > 46. These values ar e consistent with the canonical 



AGN radio-loud fraction of 10-20% (e.g.. IPonoso et al.ll2009l ). 



The 1.4 GHz radio number counts are commonly displayed as the differential counts 
dN/dS. This can be related back to the X-ray band by 

dN dN d log Lx 



dS d log Lx dS 

where 



(2) 



dN _ c Z""'™ (i$(Lx (min), 2:) df 

~(XZ. ( o ) 



dlogLx HoJ,^^^ dlogLx (l + z)2(fi^(l + z)3 + fi^)i/2 

In this last expression (i$(Lx(min), 2;)/(ilogLx is the hard X-ray luminosity function, di is 
the luminosity distance, and Lx(min) is the minimum X-ray luminosity that can produce a 
given 1.4 GHz radio flux 5*. The radio flux 5* is related to the rest-frame 1.4 GHz luminosity 
Lu by 

S = ^^d^±^ (4) 

The last factor in equation [2], dlogLx /dS, cannot be expressed analytically (especially when 
additional flux due to star-formation is included), and is calculated numerically. 

Thus, for a given S, the integral in equation [3] is calculated from Zram = to z^ax = 5 
with a new value of Lx(min) for each evaluation of the integrand. Separate calculations are 
performed for radio-quiet type 1 AGNs, radio-loud type 1 AGNs, radio-quiet type 2 AGNs, 
and radio- loud type 2 AGNs. The total 1.4 GHz dN/dS is then the sum of the results of the 
four individual calculations. 



- 8 - 



3. Results 

The solid line in Figure [2] plots the predicted Euclidean-normalized radio counts from 
AGNs at 1.4 GHz. In this figure the contributions from obscured, type 2 AGNs and un- 
obscured, type 1 AGNs are indicated by the dotted and short-dashed lines, respectively. 
Obscured AGNs dominate the predicted radio number counts from 0.1 /iJy to 10 mJy by 
factors > 3. This result is independent of the details of any of the modeling and is a very 
general conclusion. The X-ray background requires type 2 AGNs to outnumber type Is at 
all but the highest luminosities (see Figure [1]) , and this fact immediately translates to the 
radio number counts. Therefore, AGNs detected at /iJy fiux levels by the EVLA or the SKA 
will very likely be obscured type 2 objects. 

The black symbols in Figure [2] plot the observed total 1.4 GHz number counts taken 
from various published surveys over the last decade. The fiattening of the counts at fiuxes 
below 1 mJy is reproduced by every dataset. By construction the predicted counts pass 
through the majority of the data points for S > 1 mJy, as radio-loud AGNs dominate this 
regime, and the radio-loud fraction was determined by matching the model to the mJy data. 
At lower fiux densities the model passes well below the observed data, but this is to be 



expected as star-forming galaxies should dominate the /iJy radio counts (iGruppioni et al. 



20031 : iHuynh et al.ll2005l ). However, the model also significantly underpredicts (by a factor 



of ~ 5 at 30 /xJy) th e estim ate d number counts due t o AGNs (the red and blue data, taken 



from lPadovani et al.l (120091 ) and lSeymour et al.l (120081 ) . respectively). Therefore, our baseline 



model that is derived from an accurate fit to the X-ray background cannot account for the 
observed AGN number counts in the latest deep radio surveys. 

Figure [2] illustrates that the problem in the modeling arises from the contribution of 
radio-quiet AGNs. Radio-loud AGNs, denoted by the dot-dashed-line in the figure, are fixed 
by the mJy population and cannot help make up the deficit at 5" < 100 /iJy unless the radio 



l oud p opulation evolves very differently than what is currently observed. iPadovani et al. 



(120091 ) estimated the radio-quiet AGN contribution to the 1.4 GHz number counts in the 
CDF-S field, and these are plotted as the cyan data. The model radio-quiet population 
underpredicts the observed estimates by factors of ~ 3-5. 

The breakdown of the 1.4 GHz number counts into AGNs of differing X-ray luminosities 
is shown in Figure [31 When 5* > 1 mJy, the radio counts are dominated by moderately- 
luminous radio-loud AGNs with logLx > 43. The radio-quiet population at lower flux 
densities is dominated by low luminosity Seyfert galaxies. Indeed, AGNs with X-ray lumi- 
nosities logLx < 43 comprise the largest fraction of the radio-quiet number counts. This 
fact explains why obscured AGNs are the greatest contributor to the number counts at /iJy 
fluxes. 
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The dependence of the AGN 1.4 GHz counts on source redshift is shown in Figure |H 
Not surprisingly, the radio-loud AGNs are predominately found at 1 < 5; < 3 where higher 
luminosity AGNs are more common. However, once the radio-quiet population starts dom- 
inating the counts at S* < 100 /iJy, we find that Seyfert galaxies at z < 1 are the most 
common contributor to the radio-quiet AGN counts for 0.5 < 5" < 100 fiJj, but at lower 
fluxes the counts are dominated by increasingly more distant AGNs. This behavior is in 
contrast to models of the radio-counts where pushing to fai nter and faint er flux limits only 



reveals less luminous sources at roughly the same distance (jCondorull989l ). Figures [3] and H] 
show that this is not the case for radio surveys of AGNs. In fact, these calculations show 
that /iJy radio observations will be a good probe of tracing the evolution of Seyfert galaxies 
over a large fraction of cosmic time. 

The above results have shown that the predicted 1.4 GHz radio counts from radio-quiet 
AGNs are ~ 5 x too low when compared against the latest estimates from deep surveys. One 
way of increasing the model AGN radio counts is to add further Compton-thick AGNs. If the 
additional Compton-thick AGNs had very large columns with log A^h ^ 25 then, in principle, 
the X-ray background would be insensitive to their presence. However, the Compton-thick 
fraction has to be increased to 10 times the number of Com pton-thin type 2s in o rder to 



bring the radio-quiet number counts to the level estimated by iPadovani et al.l (120091 ). Even 
if all these Compton-thick AGNs has column densities logA^H > 25, the peak intensity of 
the X-ray background spectrum would be grossly overestimated. We conclude that a large 
population of hidden Compton-thick AGNs can not account for the deficit in the predicted 
1.4 GHz counts. 



3.1. The Role of Star Formation 



A natural way to increase the radio luminosity from AGNs is to have ongoing star- 
formation in the host galaxy. The fact that the AGNs contributing to the fiJj radio emission 
have significant nuclear obscuration is also a clue that star-formation m ay be important, as 
it may play a significant role in obscuring Seyfert galaxies at z ~ 1 (IFabian et al.l Il998 



Wada fc NormanI I2OO2I : [Thompson et all l2005l : iBallantvnd 120081 ). The calibration of Bell 
(I2OO3I ) is used to calculate the 1.4 GHz radio luminosity fror n a given SFR, a nd a = —0.8 is 
assumed for the radio spectrum of the star-forming regions (lYun et al.l 1200 ll ). To calculate 
the number counts with both AGN and star-formation emission, the star-forming 1.4 GHz 
flux at a given SFR and z is added to the AGN flux, and then the integral in equation [3] 
is computed as before. As type 2 AGNs dominate the AGN radio counts at all fluxes (see 
Fig. [2]), emission from star-formation was only added when calculating the total flux of 
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obscured AGNs. Calculations including star-formation in type 1 AGNs had no impact on 
the results. Likewise, although the star-forming radio flux was added to the radio-loud 
population, the AGN component easily overwhelmed the emission from star-formation and 
the counts at S* > 1 mJy did not change. 

The black curves in Figure [5] show the predicted 1.4 GHz radio counts when a constant 
SFR of 3 Mq yr~^ was assumed for all obscured AGNs, independent of Lx or z. With this 
simple addition the model is now very close to the observed estimates of both the total and 
radio- quiet AGN contributio n. This modest SFR would not be out of place in the Milky 
Way (ISolomon &: Saed Il988f). and is cons istent with constraints from the cosmic infrared 
background (IBallantyne fc PapovichI 120071 ). The fact that adding such a moderate amount 
of star-formation in AGN host galaxies can bring the models in line with the data is strong 
evidence tha t such levels of star-for mation is common in the host galaxies of obscured AGNs 
at 2 > 0.5 (ISilverman et al.ll2009l ). Moreover, the deep /iJy radio surveys will be able to 
constrain the average SFRs in obscured AGNs and determine any dependence on luminosity 
and z. 

As an illustration of what may be expected, we consider a SFR that is a function of 
both redshift and the AGN luminosity. The luminosity dependence may be expected if star- 
forming d isks at distance ~ 100 pc froni the c entral engine both obscure and feed the blac k 



hole fe.g.. lShi et al.ll2007l : IWatabe et al.ll2008f ). In the analytic models of iBallantvnd fl2008f ) 



the maximum SFR found in these disks was proportional to the X-ray luminosity, so that 
less-luminous Seyfert galaxies required a lower SFR than more luminous AGNs. The reason 
is that to fuel a luminous AGN requires a significant gas supply that can also turn rapidly 
into stars. Using these results as a guide, the luminosity dependence seems to roughly follow 
SFRoc (logLx -40)3-^ 



To estimate the redshift de pendence, we note that iKim et al.l (120061 ) mention that the 
SDSS type 2 AGNs identified bv lZakamska et all J2003h have a SFR~ 20 Mq yr'^ as judged 
by the [O II] lii ie strength. The me an [O III] luminosity for their sample is L[o m] ~ 



3 X 10^^ erg s ^. iHeckman et al.l (120051 ) found that for local type 1 AGNs, the average ratio 



between Lx and L[o m] is 1.59 dex. These authors do not correct for absorption in their 
sample of type 2 AGNs to find the ratio for ob scured AGNs, so we fol low the AGN unified 
model and apply the type 1 conversion to the Zakamska et al. ( 2003h sample resulting i n 



Lx ~ 10^^ e rg s \ The majority of the SDSS objects are at 2; ~ 0.3 (IZakamska et al.ll2003l ). 



Ballantynd (120081 ) found a SFR~ 40 Mq yr~^ in a nuclear starburst disk fueling an AGN 
with Lx ~ 10^"^ erg s~^. If we assume that this situation is common at z ~ 1 (i.e., at the 
peak of the cosmic star-formation history), then the SFR in the AGN host galaxy varies as 
(1 + zY''^^. Surprisingly, this is comparable to the redshift dependence of the SFR in quasar 
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host galaxies inferred by ISerjeant fc Hatziminaogloul (120091 ) using a sample selected in the 



far-infrared; therefore, this redshift dependence does not seem unreasonable. As with the 
AGN type 2/type 1 ratio, the redshift evolution was halted at z = 1. 

The normalization of the SFR evolution was adjusted to obtain an adequate looking fit 
to the radio data. The final evolutionary law is 

SFR ^ 0.25(1 + 2)^-^*^(logLx - 40)^-^ Mq yr-\ (5) 

for z < 1 (the SFR remains at the z = 1 values at larger redshifts) and is plotted in Figure [61 
The redshift and luminosity dependence s in equation O are steeper than those that were 



recently found by [Silverman et al.l (120091 ) by studying AGN hosts in the zCOSMOS survey. 
However, the scatter in the data is about an order of magnitude, and the predicted SFRs 
fall within the range observed in the zCOSMOS data. Given the uncertainties in both the 
modeling and observations, equation [5] should be considered as a reasonable, observationally- 
motivated suggestion for how the SFR in AGN host galaxies may evolve. 

The 1.4 GHz radio counts from AGNs and their host galaxies, assuming a SFR following 
eq.[5l are shown as the green curves in Figure[5l The new predicted radio-counts pass through 
the majority of the observed data points. AGNs with logLx < 43 still provide the majority 
of the counts at < 100 /iJy, but higher luminosity AGNs are now a significant contributor 
at S" ~ 40 /xJy (Figure [7]). Comparing Figs. [7] and [21 we see that the star- formation has the 
largest impact on the counts only for AGNs with log Lx < 44. This is because the radio-flux 
from the AGN grows faster with luminosity than the flux from the star-forming regions. The 
/iJy counts will then be most sensitive to measuring the SFR in Seyfert galaxies. 

With the addition of the redshift-dependent star-formation law, the redshift distribution 
of the AGN radio-counts spreads to higher fluxes (see Figure [HI). In this model, radio-quiet 
AGNs at z < 1 are predicted to be common at S" ~ 40 /iJy, as compared to the case with 
no star- formation where these AGNs peak at ~ 10 fiJj (Fig. [1]). The planned EVLA 
surveys should measure the total radio-counts down to ~ 1 /iJy. Therefore, the evolution of 
star- formation in AGN host galaxies should, in principle, be measured by the deep EVLA 
surveys. 



3.2. Is Star Formation Required? 

The results presented above indicate that additional radio emission from star-forming 
regions in AGN host galaxies must be included in order for the predicted number counts to 
match the observations. However, this conclusion will depend on the choice of the X-ray to 
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radio conversion employed in the modeling. We have argued i n Sect. [3] that a luminosity- 
dependent conversion factor (eq. [p drawn from t he wo rk of iTerashima &: WilsonI (120031 ) 
was the most appropriate choice. IPadovani et al. ( 2009 ) found that their estimates of the 
radio-quiet contribution to the number counts were adequatel y described by a model w ith 
a luminosity-independent conversion factor consistent with the iBrinkmann et al.l (120001 ) re- 
lation. To test the sensitivity of ou r conclusions we the refore computed a new model of 
the 1.4 GHz radio counts using the iPanessa et al.l (120071 ) relationship derived from high- 
resolution observations of local low-luminosity Seyfert galaxies: logLx ~ (0.97) log z/Li,(5 GHz) + 
5.23. Figure [9] presents the results of this calculation and clearly shows that this model does a 
good job fitting the estimated AGN counts with no additional emissio n from star-formation . 
The contribution from radio-quiet AGNs is slightly overpredicted, but lPadovani et al.l (120091 ) 
state that size of the error-bars are most likely underestimated. While emission from star- 
formation in t he AG N host galaxies is not required to match the data while using the 
Panessa et al.l (120071 ) relation, a small (~ 2 Mq yr~^) amount of star-formation in the host 
galaxies is still consistent with the observational constraints. 

Figures fTOl and fTTl plot the contributions to the counts from AGNs in different luminosity 
and z bins in this scenario. A comparison between these figures and Figs. [7] an d [5] show signif- 



i cant differences in the properties of AGNs that dominate the /xJy counts. The lPanessa et al. 
(120071 ) conversion results in larger radio luminositie s for AGNs with log Ly ~ 43-44 and, as 



these AGNs dominate the X-ray background (e.g.. iBallantyne et al.ll2006l ). they also domi- 
nate the yuJy radio counts. This is in contrast with the model that requires star-formation 
where AGNs with logL^ ^ 43 will be most common at these fluxes (Fig. [7]). Similarly, 
Fig. [11] predicts that AGNs at 1 < 2; < 3 will be equally common at /iJy fluxes as AGNs 
at ^ < 1 , while Fig. [8] shows that AGNs at z < 1 will dominate if there is significant star- 
formation within the host galaxies. Therefore, the luminosity and redshift distributions of 
the AGNs detected in the deep radio surveys may help indicate whether or not additional 
radio emiss i on fro m star-formation is necessary to explain the count distribution. Recently, 
Tozzi et al.l (120091 ) has examined the properties of X-ray sources in the C handra Deep Field 



South that were a l so det ected by the deep VLA survey of that region (iKellermann et al. 



20081 ) . iTozzi et al.l (120091 ) found that the redshift distribution of the radio sources with X- 
ray counterparts was peaked at lower redshifts (median 0.73) than the X-ray sources without 
a radio match (median 1.03). This result is qualitatively consistent with the predictions of 
the AGN number counts that includes star- formation from the host galaxy (Fig. [8]). 
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4. Discussion and Summary 



The principle conclusion of this work is that, depending on the nature of the X-ray to 
radio conversion in radio-quiet AGNs, the observed radio number counts of AGNs indicate 
that moderate levels of star-formation exist in the radio-quiet fiJy AGN population. Thus, 
deep radio surveys will be an important probe of the connection between galaxy and black- 
hole growth. In particular, combining /iJy radio data with results from surveys at other 
wavelengths will probe how the star-formation properties of AGN host galaxies evolve with 
z. 

The model of the radio-counts was based on a fit to the X-ray background spectrum 
and made use of the hard X-ray luminosity function. It was found that in almost all cases 
the /iJy counts due to AGNs are dominated by Seyfert galaxies with logLx < 43. The radio 
flux produced by more luminous AGNs was swamped by the minority radio-loud population. 
The implication is that AGNs selected by their fiJj radio emission will be overwhelmingly 
obscured Seyfert galaxies spread over a range of redshifts. Radio selection, therefore, will 
result in a way of probing the evolution of the low-luminosity end of the hard X-ray lumi- 
nosity function. This will be especially interesting in light of recent suggestions that the 
evolutionary pathways may be ve ry different for galaxies that host high-luminosity quasar s 
and those that host Seyferts (e.g.. lBallantyne et al.ll2006l : lHasingeiil2008l : iHickox et al.ll2009l ). 
The idea here is that quasars are the signatures of the formati on of massive galaxies gener 



Kaufmann fc Haehnelt 
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ated by the violent m ergers of two large gas-rich galaxies (e.g 
Hopkins et al.ll2005l ). The lower-luminosity AGNs, such as Seyferts, do not require such 
large fueling rates, and may be signatures of a less violent path of galaxy assembly. Thus, 
the Seyfert galaxies found at /iJy radio fluxes will be a simple way to select and measure 
the properties of AGNs and their host galaxies with little contamination from the higher 
luminosity population. 

The results described in Sect. [3] show that moderate amounts of star- formation may 
be ongoing in the /iJy AGN population. In addition, these AGNs will be overwhelmingly 
obscured, or type 2 AGNs. This is consistent with the i dea that the abso rbing medium 
around an AGN is related to the host galaxy SFR (e.g., iBallantynd l2008l ) . A corollary 



of this scenario would be that the SFR and its associated radio-fiux would be correlated 
with X-ray column dens i ty. T here are some observ a tional hints that this may be occurring: 
both iGeorgakakis et al.l (120041 ) and [Richards et al.l (120071 ) found an increase in the number 
of X-ray sources with a radio counte rpart in a deep ra dio survey as the X-ray AGNs became 
progressively more obscured (but see lTozzi et al.ll2009l ). Breaking down the /xJy radio counts 
into contributions from different X-ray obscuring columns requires a detailed model of the 
circumnuclear star-forming region and will be the subject of future work. 
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In summary, the increased radio continuum sensitivity of the EVLA and the SKA will 
prove extremely powerful in the study of AGN evolution. Models of the 1.4 GHz number 
counts derived from a fit to the hard X-ray background show that the /iJy population of 
AGNs will be dominated by radio-quiet, obscured Seyfert galaxies at 2; ~ 1. In addition, 
the radio fluxes from these AGNs may be enhanced by galactic star-formation. Indeed, 
current estimates of the AGN contribution to the 1.4 GHz radio counts seems to require 
moderate levels of star-formation in order for the models to match the observed data points. 
Multiwavelength foUowup observations of AGNs selected by their //Jy radio emission will be 
able to trace the evolution of Seyfert galaxies from 2; ~ 3 to the present day. 

The author thanks V. Smolcic, P. Padovani, N. Seymour and F. Bauer for useful con- 
versations and comments. The anonymous referee is acknowledged for a detailed and con- 
structive report. 
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Fig. 1. — The fraction o f type 2 AGNs as a function of Lx as determined by the Swift /BAT 
survey of iTueller et al. ( 2008 ). The observed correlation between Lx and -^14-195 kev 
(jWinter et aLl 120091 ) was used to convert the luminosities to the Lx band. The solid line is 
a power-law fit to the data assuming a 4:1 type 2/type 1 ratio at logLx = 41.5. The best 
fit has a reduced = 0-97 and follows /2 oc logL^^^"'^. 



- 21 - 



10^ 



10 - 



1 - 



0.1 ^ 



§ 10-2 
in 



10- 



10- 



VLA-COSMOS 

* CDF-S (AGN.RQ AGN 
A SSA13 

X la*"- XMM-Newton (AGN) 

* Suburu/XMM-Newton Deep Field 

* VLA-VVDS 
A Phoenix 

□ FIRST 




--- Type 2 AGN 
- Type 1 AGN 

— Radio Quiet AGN 4 

- Radio Loud AGN 



0.1 



1 10 100 1000 

1.4 GHz Flux Density (/^Jy) 



10^ 



Fig. 2. — The solid line plots the predicted Euclidean- normalized 1.4 GHz radio counts 
from 0.1 /iJy to 10 mJy obtained from the X-ray background modeling described in Sec- 
tion [21 The dotted and short-dashed lines plot the contribution to the counts from ob- 
scured and unobscured AGNs, respectively. These curves include AGNs of all radio pow- 
ers. Type 2 AGNs dominate both the radio-loud and radio-quiet populations at these 
flux levels. The contribution from radio-quiet AGNs is plotted as the long-dashed line, 
while radio-loud AGNs are shown as the dot-dashed line. These curves include both ob- 
scured and unobscured AGNs. The points indica te the obs erved counts obtained f rom var- 
i ous surveys: V. 



(IFomalont et al 



20061). 13hr XM M-Newton fISevmour et a 



A-COSMOS fiBondi et al.ll2008D. CDF-S flKellermann et al.ll2008[ ). SSA13 



20031) and FIRST (IWhite et al. 



Deep Field flSimpson et all 12006 ). VLA-VVDS flBondi et al 



20041) . Suhmu / XMM-Newton 



20031 ). Phoenix (IHopkins et al. 



19971 ). The blue and red data plot the estiniated AGN 



contribu tion to the radio cou nts from the 13hr XMM-Newton (jSeymour et al.l 120081 ) and 
CDF-S ( jPadovani et al.ll2009l ) surveys, respective ly. The cyan data po ints are estimates of 
the radio-quiet AGN contribution in the CDF-S (jPadovani et al.l 120091 ). 
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Fig. 3. — As in Figure [2l but the contributions to the total number counts from AGNs with 
different absorption-corrected X-ray luminosities are now indicated. 
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Fig. 4. — As in Figure Ej but the contributions to the total number counts from AGNs at 
different redshifts are now indicated. 
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Fig. 5. — The black solid line plots the predicted number counts from AGNs including radio 
emission from a constant SFR of 3 Mq yr~^. The black dashed line shows the contribution 
from radio-quiet AGNs in this scenario. The green lines are the analogous curves for a model 
where the SFR is a function of both logLx and z (eq. [5]). In both cases the radio emission 
was added to only the type 2 AGNs. The data points are the same as in Fig. [21 
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Fig. 6. — Contours of SFR in Mq yr~^ as a function of logLx and z calculated using 
equation [5l This law indicates that SFRs < 30 Mq yr~^ may be common in the host 
galaxies of the /xJy AGN population. 
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Fig. 7. — As in Figure [3], but radio flux from a SFR given by eq. [5] was added to tlie type 2 
AGNs. 
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Fig. 8. — As in Figure HI but radio flux from a SFR given by eq. [5] was added to tlie type 2 
AGNs. 
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Fig. 9. — The solid l ine plots the predicted number counts due to AGNs when assuming the 
Panessa et al.l (120071 ) relationship between Lx and vLyi^ GHz). No additional emission from 
star-forming regions has been included. The contributions from radio-quiet and radio-loud 
AGNs are shown as indicated. 
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Fig. 10. — As in Figure [9|, but the contributions to the total number counts from AGNs with 
different absorption-corrected X-ray luminosities are now indicated. 
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Fig. 11. — As in Figure [9l but the contributions to the total number counts from AGNs at 
different redshifts are now indicated. 



